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ABSTRACT 



The elastic transient response of an impulsively loaded circular 
ring is examined. The numerical stability criterion associated with 
a finite difference formulation of the problem is briefly discussed, 
and a method to reduce the stability restriction is presented. This 
method assumes that the density of the ring is greater in the tangen- 
tial direction than it is in the radial direction. The analytical 
solution to the ring equations with the artificial tangential inertia 
is determined for a specific impulsive load. This solution is 
compared with the solution to the ring equations with the correct 
tangential inertia. The comparison shows that on the basis of pre- 
dicting peak values of stress, the results for the ring with the 
artificial tangential inertia are very close to those of the unmod- 
ified ring. Thus, the method reduces the numerical stability require- 
ments and does not significantly affect the principal results. 
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CHAPTER I 



INTRODUCTION 

The finite difference method of analysis has long been 
recognized as an effective and practical means for solving partial 
differential equations on a digital computer. The accuracy asso- 
ciated with this method can be adjusted by changing the incremental 
sizes of the independent variables involved. In the limit, when 
the sizes of the increments approach zero, the solution to the 
finite difference equations becomes identical to the solution to 
the corresponding differential equations. Consequently, the engineer 
may improve the accuracy of his finite difference solution by choos- 
ing sufficiently small values for the incremental sizes of the inde- 
pendent variables, with greater accuracy accompanied by increased 
execution time on the digital computer. However, there are many 
problems where the finite difference method imposes a restriction 
on the sizes of the increments. This restriction is known as the 
stability criterion and arises from a purely numerical phenomenon. 

The stability criterion defines a critical limit for the ratio of 
the increment sizes above which a numerically stable, valid solu- 
tion will result. A value of the ratio below this limit will lead 
to a divergent, meaningless solution. 

The stability criterion for the finite difference dynamic 
analysis of the thin, homogeneous, elastic, circular ring has been 
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It can be simplified to 



defined by Ball. ^ 

AA4> . . D ' , , 

AT "y ("h ( I 

where A is the radius of the ring; h is the thickness of the ring; 

D = Eh, the stretch rigidity; ^ is the ring density; and, A 
and AT are the angular and time increments, respectively. A 
meaningful numerical solution for the transient response of the ring 
is possible only when the above restriction on A<p and AT 
is observed. A physical interpretation of the stability criterion 
says that the finite difference solution must propagate around the 
ring at least as rapidly as a stress wave. 

Once the size of a 4 > is chosen, the maximum size for the 
time increment becomes 

= (a4>) 

An accurate solution for the ring response requires a small value 
for the angular increment, . Therefore, according to Equation 

(1), ( AT is automatically restricted to a small value; and 

hence, large computation times will be required on the computer. 

To reduce the restrictions on the numerical stability of the 
ring response a method is proposed whereby the maximum time increment 
may be increased without producing an unstable solution. According 
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to Equation (1), ( AT is determined once is chosen. 

However, if the density ^ is considered to be a variable in the 
stability criterion then ( AT will be proportional to the 

square root of the density at a particular value for A ^ • The 
density in Equation (1) is a measure of the ring’s inertia in the 
direction of the propagation of stress disturbances around the ring, 
that is, the tangential direction. The method proposed here modifies 
the ring by introducing an artificial tangential density that is 
larger than the actual ring density. Thus, the ring inertia will be 
characterized by the density ^ in the radial direction and by a 
different density, ^ ^ , in the tangential direction. At a given 

value of A <p y the maximum time increment for a stable response 
of the modified ring will be proportional to the square root of 
and, hence, will be larger than the critical time increment for the 
original ring by the factor /r. 

For example, for ^ 4 the time increment for stable response 

of the modified ring is twice as large as the time increment for the 
original ring for the same A(^ • Thus, the computation time can 

be cut in half. The decreased computational time afforded by the 
proposed method is a significant improvement on the conventional 
finite difference method of analysis. However, the introduction 
of the artificial tangential inertia will also lead to an incorrect 
solution for the ring response. Thus, the value of the method will 
depend upon the accuracy of the solution of the modified ring. 
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This paper investigates the effect of an artificial tangential 
inertia on the analytical solution for a particular ring problem. 
Solutions are obtained for the elastic response of an impulsively 
loaded circular ring with and without artificial tangential inertia. 
The maximum stress in the ring is computed for the original ring 
( y * 1) and for rings with * 2, 3, and 4. A comparison of the 
maximum stress for each ring demonstrates the value of the method. , 
This study was completed during the 1967-1968 academic year at 
the Naval Postgraduate School, Monterey, California. Acknowledge- 
ment is gratefully made to Professor Robert E. Ball of the Aeronautics 
Department for his guidance and consultation. 
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CHAPTER II 



A DISCUSSION OF THE RING EQUATIONS 



This thesis considers a thin, circular ring composed of a 
homogeneous material. The governing equations are derived in 
Reference 4 and are repeated here for completeness. Figures I and 
2 illustrate the general geometrical characteristics. The radius A 
is measured to the middle surface of the ring, which lies at a dis- 
tance — from both edges. A differential element on the ring, d , 

is located at an angle ^ from the bottom of the ring, where is 

measured positive counter-clockwise. The displacement of a point 
on the middle surface is described by a radial and a tangential 
component, W and V, respectively. The radial distance from the 
middle surface to some arbitrary point on the ring is given by Z in 
Figure 2. 

A two-dimensional problem will be assumed with regard to the 
internal stresses of the ring. Figure 3 illustrates the forces 
acting on the element d4> due to the normal and shear stress dis- 
tributions on the cross-sections. The forces N, M, and Q are, 
respectively, the normal force, the bending moment, and the radial 
shear force. These forces are stress-resultants referred to the 
middle surface of the element. The definitions of N and M are given 
by the following expressions 





N 



cr^ dz 



£ 



(2a) 




J.A 



2 
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FIG. I 

GENERAL RING GEOMETRY 
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FIG. 2 

RING MIDDLE SURFACE 




FIG. 3 

TWO-DIMENSIONAL STRESS-RESULTANTS 
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M = 



h 

z 



li 






z dz 



/-ii 

2 



zdz 



Ul 

2 



(2b) 



where and are the stress and strain on the cross section 

at a distance Z from the middle surface. 

In order to obtain an expression for the strain in terms 

of Z, two additional assumptions are made: 

1. A line normal to the middle surface before deformation 
will be normal to this surface after deformation. 

2. A point on the normal to the middle surface will have a 
constant distance Z from this surface, both before and 
after deformation. 

Accordingly , 



< - / J #w w 

^ 'A A^z 

Substituting Equation (3) into the definitions for N and 
M, Equations (2a) and (2b) , leads to 




^<P 






M 



K I aV 
AM 



■h w 



(4) 

(5) 
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where 



£/? 



- K ^2 

The linear equilibrium equations for the element follow directly 
from an examination of Figure 3. They are 



ML 

^<P 



-Q - 



+ /V = |oAh 






(6a) 



(6b) 



= 0 



A ^<t> 

where rotary inertia is neglected. 

Eliminating the radial shear Q from these equations leads to 



(6c) 



M. _ _L_ JM. 



A ^<p 



I 



iM 



A 



i + A/ - 



f>hh 






(7a) 

(7b) 



Substituting Equations (4) and (5), relating the s tres&resul tants 
to the displacements W and V, into Equations (7a) and (7b), gives 



_D_/i!v . AW 

A ( 

+ iv] 






AM ^<#i" 






iV 



(8a) 

(8b) 



These are the differential equations that describe the free 
vibrations of the thin, two-dimensional, elastic ring. 
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CHAPTER III 



SOLUTION FOR THE IMPULSIVELY LOADED RING 

In order to adequately test the proposed method for increasing 
the critical time step, an impulsive load that causes both large 
bending and membrane stresses should be employed. Consequently, 
the impulsive load illustrated in Figure 4 is taken as the loading 
on the ring at time t = 0. The impulse is applied to the top half 
of fhe ring and is characterized by a cosine distribution directed 
radially inward. 

The initial conditions for the radial cosine impulse are given 
by 



w , 0 ) = 0 v( 4 >>o) = 0 

^{<*>, 0 ) =0 

-|^ (4>.0) = -V, C0S<(> 

= 0 (9) 

where Vo = 

Since the ring is unrestrained, a net momentum will be imparted 
by the radial velocity acting on the upper half of the ring. This 
will cause a rigid body translation of the ring, and hence, unbounded 
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FIG. 4 

IMPULSE LOADING ON THE RING 
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values for the displacements V and W. To avoid this situation, the 
initial conditions must be modified in such a way that the rigid 
body motion is eliminated, while preserving the internal effect of 
the impulse. Reference 2 defines the following modified set of 
initial conditions for this loading 



W(<^; 0)=0 



V'(^; 0} =0 



^ SIN 

^(4,, 0) = -;|-l4 cos 

^ ^ cos <h 

4 ^ 



TT , . , 37T 



2L 

Z ^ 2. 



( 10 ) 



with this set of initial conditions, the ring will remain fixed in 
space; and hence, the displacements V and W will measure actual 
deformations of the ring. 

The following solutions for W and V are assumed 



W(4>. t) = ) A„(7) Cos n4> 



n= o 



OO 

V {<P>T) B„(r) SIN n<t> 



n=i 



A„(T)= 






i 



(lla) 



(11b) 



(12) 



*a and b are arbitrary constants, 
n n 
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Based upon these solutions and the given initial conditions the 
expressions for the bending moment and normal force are* 



and are non-dimensional representations of the bending moment 
and the normal force, both of which may be solved without specifying 
the impulse magnitude, . The other terms will be described later • 
Equations (13a) and(13b) characterize the dynamic response of the 
elastic ring due to the specified impulse. Before proceeding further, 
however, a physical interpretation will be given to these expressions. 

*These expressions are derived in Appendix A. Since this is 
the original ring, J = 1. 





(13a) 



^ axT + CO^ 4> 




•h S/W 3n o>or COS Mo 



(13b) 



where 
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Consider, first, the various mode shapes excited by the impulsive 
loading. Figure 5 Illustrates the lower modes of vibration for the 
ring. The lowest mode, n 0, is characterized by a radial displace- 
ment only, where W is the same for all points. A small change in 
curvature in addition to the middle surface straining will result; 
and hence, and will both have contributions from this mode as 
can be seen from Equations (4) and (5). The next mode, n * 1, 
involves middle surface straining only. The shape of the ring 
remains the same, with no change in the effective curvature of the 
ring. Hence, this mode will contribute only to the normal force, 

No • For n > 1, only the even mode shapes appear in the expressions 
for the s tress-resultants I This situation arises from the modified 
initial conditions imposed on the ring. At time t ® 0, positive 
radial velocities are required at the top and bottom of the ring, 

^ 7T and 0, while zero radial velocities are prescribed at the 
sides, 4^* ^ * Referring to Figure 5, it can be seen that only 

the even modes can account for these velocities and hence, are 
the only ones excited by the Impulse. 

Finally, examination of Equations (13a) and (13b) reveals th^t 
two frequencies, ^ associated with each mode 

of vibration of the ring. These frequencies arise from the fact that 
the ring element, , may displace independently in two separate 

directions. From Appendix A, the frequency coefficients, and 
g^, are approximately given by 
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FIG. 5 

LOWER MODE SHAPES 
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fy, - n(n"-/) 



(lA) 



3n = (rFhT 

Since a thin ring is assumed, is very small, and f << g . The 

n n 

frequency coefficients, f^ and g^, characterize, respectively, the 

radial and tangential vibrations of the ring. Due to the relative 

magnitudes of f^ and g^, the expressions for and will each be 

dominated by one of the two sine terms of the series. Hence, 

will vary at a frequency given approximately by ^ while 

will be characterized by a much larger frequency % 

To obtain numerical data from Equations (13a) and (13b), it is 

necessary to select physical data describing the ring. Figure 6 

outlines the physical characteristics for a thin ring composed of 

aluminum. The frequency coefficients, f and g , have been calcu- 

n n 

lated on a digital computer at 20 micro-second intervals over a 
period of time corresponding to one complete cycle for . Stresses 
on the inner and outer edges of the ring were also calculated from 
the following non-dimensional relationships 

(outer edge) (jJ * 

t Vo 

(inner edge) O’.' = 

t. V© 



= N, ^ M, (15a) 

= N, - tSt Me (15b) 
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16980 rod. /sec. 


K = 


1.43 X 


10^ Ibf.-in. 


2.74 X 10* Ibf./in. 




3.14 


X Id* 

m. ^ 


10.9 X 10* Ibf. /in.* 


a « 


3.54 


X 10'® 



FIG. 6 

PHYSICAL RING DATA 
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The maximum values for , and CT^ found over the time sampled 

have been listed in Table 

Figures 7 through 10 illustrate the dynamic response of the 
ring through the first half cycle of It should be noted from 

Table I that the response at cp s: ^ 'SiL is of secondary importance 

and has been omitted from the curves. Also, the response at 4^= 7T 
is identical in character to that at <^ = 0 and has been omitted. 

The symmetry of the even modes make it unnecessary to determine 
solutions for ring locations in the range ^ <p ^ 27 T. 



*The maximum was observed to be within a 

the maximum orj found. 



few per cent of 



positive indicates a negative bending moment. 
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TABLE I 



MAXIMUM VALUES FOR THE STRESS-RESULTANTS 















M 


G-o 

MAX 



0 


2.22 


-129.6 


-3.45 


IT 


1.80 


26.5 


1.99 


4 








7T 


1.31 


125.8 


2.57 


2 








37T 


1.77 


26.5 


1.80 


4 








7T 


2.39 


-129.6 


-3.27 
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I40^ 
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BENDING MOMENT RESPONSE AT ^ = 0 (y=l) 
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FIG. 8 

NORMAL FORCE AND EDGE STRESS AT = 0 (y =1) 



14 0 
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FIG. 9 

BENDING MOMENT RESPONSE AT <f> = 
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NORMAL FORCE AND EDGE STRESS AT 



CHAPTER IV 



SOLUTIONS FOR THE IMPULSIVELY LOADED 



MODIFIED RING 



The assumption that the ring’s inertia characteristics are 
variable with regard to tangential vibrations will now be considered. 
The ring equations must be modified to account for this assumption, 
and modified solutions analagous to the original ring are obtained 
in Appendix A for the impulse loading described in Chapter II. 

Figure 11 illustrates the assumed inertia forces associated 
with the ring element, d Cp • A simple adjustment of the ring 
equations will account for this situation. The following expres- 
sions for the stress -resultants , derived in Appendix A, are obtained 
for the modified ring 



Mo = SIN cOot -h 







-h SiN COS ncp Mo 
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FIG. II 

ASSUMED INERTIA FORCES 
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consider 



To determine the effect of the factor ^ on and ^ , 
first the change in the frequency coefficients, f^ and As deter- 

mined in Appendix A, these coefficients are 



-Pn = n(n"-/) 



o< 






( 18 ) 



9n - 






( 19 ) 



For this analysis, it will be assumed that takes on values between 
1 and 4 * Therefore, f^ will not be significantly affected by 6 , 
especially for the higher mode shapes. However, g"^ will be very 
sensitive to a change in 5 , as can be seen in Equation ( 19 ). Figures 
12 and 13 graphically illustrate the error in the modified frequency 
coefficients as a function of . A physical interpretation of the 
factor y will verify the above results. As ^ increases, a progres- 
sively heavier system is seen in the tangential direction, giving 
rise to a response of lower frequency, • The radial vibrations, 

characterized by f , should not be significantly altered by 

n ^ 

since the density in this direction is always given by the original 
ring density, ^ . 

Values for the modified frequency coefficients have been calcu- 
lated and listed in Appendix B for the ring defined in Figure 6. 

Again, it can be seen that ^ K > then, will be 

approximately described by sinusoidals at frequencies > 

respectively. 
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FIG. 12 

CHANGE INT„ with 7 



1.0 



9n 

9n 





' I I 

2 4 6 8 n 

FIG. 13 

CHANGE IN WITH 7 
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Examination of the expressions for and reveals the 
dominant terms. For the bending moment, Equation (16), the lower 
frequency sine term will dominate. The magnitude of the coeffi- 
cient for this term will decrease as the product, ^ g increases, 
corresponding to increasing values for • Thus, the time 
response for the bending moment will decrease in magnitude as 
increases. The wave form for the response, however, will not change 
significantly. Examining Equation (17) in a similar manner, it can 
be seen that both the wave form and the magnitude of the time response 
for No will change as ^ increases. 

Figures 14 through 19 illustrate the curves for , No , and 
0^^^ for three values of at <p ^ 0 • Figures 20 through 22 show 
Kq for ^ values obtained for and N*o at 4^ * ^ ^-tid 

^ were found again to be small compared with other locations 

on the ring and have been omitted; Also, the curves for * TT 
have not been included because of the similarity with the response 
for cp> ^ 0, 

Tables II through IV list the maximum values observed for the 
solutions of and (tJ . The corresponding values for the 

unmodified ring have also been included for comparison. 



and crj at 

magnitudes. ^ 



are not Included because of their small 
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TABLE II 



MAXIMUM BENDING MOMENTS 



4 > 


^=1 


jr=2 


y=3 


^=4 


0 


-129.6 


-129.5 


-124.9 


-117.7 


7T 

2 


125.8 


118.3 


110.4 


103.2 


7T 


-129.6 


-129.5 


-124.9 


-117.7 






TABLE III 






MAXIMUM 


NORMAL 


FORCE (MAGNITUDES) 


4 > 




y=2 


ar=3 




0 


2.22 


2.38 


2.50 


2.48 


7T 

2 


1.31 


1.48 


1.14 


1.13 


7T 


2.39 


2.43 


2.46 


2.48 





TABLE 
MAXIMUM OUTER 


IV 

EDGE STRESS 




JT- / 


jr=2 






0 


-3.45 


-3.49 


-3.28 


3.57 


IT 

2 


2.57 


2.38 


2.20 


2.18 


7T 


-3.27 


-3.09 


-3.31 


2.98 
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I40^ 
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FIG. 14 

BENDING MOMENT RESPONSE AT«^=0(y = 2) 
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FIG. 15 

NORMAL FORCE AND EDGE STRESS AT = 0 (/=2) 



I40^ 
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FIG. 16 

BENDING MOMENT RESPONSE AT<^=0(/=3) 
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FIG. 17 

NORMAL FORCE AND EDGE STRESS AT </> =0 (/= 3) 



I40H 




BENDING MOMENT RESPONSE AT ^ = 0(y*4) 
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FIG. 19 

NORMAL FORCE AND EDGE STRESS AT <^»0 (/» 4) 



I40^ 
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FIG. 20 

BENDING MOMENT RESPONSE AT«^=-“ (y=2) 
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BENDING MOMENT RESPONSE AT <^ = y (y =3) 



140 -^ 




48 



FIG 22 

BENDING MOMENT RESPONSE AT f =f- (y 



CHAPTER V 



CONCLUSIONS 

A comparison of the response of the original ring with the 
results for the modified ring with ^ = 2, 3, and 4 can be made using 
the data presented in Tables I through IV and the response curves 
shown in Figures 14 through 22, A study of the time response for 
the bending moment, and Mq > reveals that the solutions to the 
modified ring, in general, range from fair ( ^ - 2) to poor ( S' = 4) 
approximations to the solution for the original ring. Average errors 
range from 10 per cent for S' = 2 , to about 40 per cent, for 8^ = 4. 
However, the percentage error associated with each modified curve 
is nearly constant, and the wave form is identical to the unmodified 
curve. It would appear, then, that the data from a modified solu- 
tion could be increased by some constant factor to yield a very 
good approximation over the first quarter cycle. 

Based upon the above comparison, the proposed method appears to 
introduce unacceptable errors into the solution for J > 2. However, 
the actual time response for the bending moment is of little signifi- 
cance in a dynamic analysis, since the maximum value of bending 
moment is usually the desired result. The maximvim bending moments 
observed for all solutions have been listed in Table II. Based on 
these values the solutions for the modified ring are very good 
approximations to the results for the unmodified ring. Errors in the 
maximum bending moment of about 4 per cent and 8 per cent can be 

seen for S' values of 3 and 4, with essentially no error observed 
at y = 2. 
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The effect of a heavier system in the tangential direction 
is readily observed in an examination of the various normal force 
response curves. Both the wave forms and the maximum values are 
affected by a change in the factor if . As anticipated, the tangen- 
tial vibrations decrease in frequency as H increases. This can be 
seen from the spacing of the beats. Despite the obvious changes 
in the normal force curves for the modified ring, the maximum peak 
values, listed in Table III, are not significantly altered as ^ 
changes . 

The primary concern of any dynamic analysis is usually the 
maximum stress that occurs in the ring. Since both the bending 
moment and the normal force contribute to the stresses, the 
maximum value will be found on one edge of the ring. The relative 
importance of each contributor, and , can readily be seen 
from the stress curves found in Figures 14 through 19 . Table IV 
lists the maximum non-dimensional values found from these curves. 
The bottom of the ring, ^ = O , is seen to be the critical stress 
location. A non-dimensional value of - 3.45 is the maximum stress 
calculated from the solution for the original ring. In comparison, 
the modified ring solutions have values of - 3 . 49 , - 3 . 28 , and - 3 . 57 , 
all of which are within 5 per cent of the original ring result. 
Consequently, the assumption of an artificial tangential density 
has a negligible effect on the ring’s maximum dynamic response for 
the impulsive loading condition considered. 



50 



In conclusion, the proposed method for reducing the stability 
restrictions presented in Chapter I appears to be valid. The same 
method when applied to the finite difference approach should produce 
results very similar to those presented here. Thus, because of the 
assumed artificial tangential inertia, an accurate finite difference 
solution should be obtained with the computation time reduced by 
about 30 per cent for ^ = 2, 40 per cent for 3, and 50 per cent 
for 4 from that of the unmodified ring analysis. 
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APPENDIX A 



SOLUTION FOR THE SPECIFIC IMPULSE LOADING 
The derivation for the solution to the radial cosine impulse 
is outlined in this appendix. It should be noted that the develop- 
ment is identical to that of the exact solution when the factor 

y = 1. 

Assuming that the density in the tangential direction is larger 
by a factor than that associated with the radial direction, the 
following modification in the ring equations will result 



A 

A’ 









(A-l) 



(A-2) 



Solutions for W and V were assumed earlier in Equations (11a) 
and (11b) and are repeated below. 



w 



oo 

(4>>7) A„(t) cos n4> 



K) = o 
oo 



v(4= ,t)= Y b„(t) siii n4> 






(A- 3) 






The symmetric and anti-symmetric forms for W and V arise from the 
symmetry of the impulse and the definition of the V, W coordinate 
system. 
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Substituting the assumed solutions into the modified ring 
equations, Equations (A-1) and (A-2) , and setting the two algebraic 



expressions equal to zero for each value of n lead to 

D(A„-r)6„) + =o (a-s) 

for n " 1 , 2 , . . . , oo 

Harmonic forma are assumed for A and B 

n n 

^ (A-6) 

B„(T) - ] 

where a and b are arbitrary constants. Substitution of these 
n n 

assumed solutions for the time functions into Equations (A-4) and 

(A-5) yields the following two expressions from whlcli the common 

i Ci>M T 

time term, e ” , has been eliminated 

On ( cly, -mhn) -h ^ A^f>h ^ 0 <a-7) 

a^, ^ o (a-8) 

The coefficients a and b cannot be determined from the above 

n n 

equations since a third unknown, COff > 1^ involved , However, they 
may be expressed in terms of one another and eliminated from both 



expressions. 



From Equation (A-7) 



= bi 



n 



n - 



Cdn 

Dn 



(A-9) 



Replacing a in Equation (A-8) with the expression above, the 
n 

constant b is common to all terms and may be factored out. A 
n 



single equation in the unknown 60^ results 



X 

y) 












j_ ^ yo< 

Y\ n 



+ n 






00c 



„ 2 

4- (v) -/) = 0 (A-10) 



where 








o< = 



K 



The modal frequency 60 for the time functions A and B may 

n n n 

now be determined. Although four values for 60^ will satisfy 
Equation (A-10), only two, the positive roots, have physical signifi- 
cance. These roots indicate that each mode of vibration is character- 
ized by two distinct frequencies 



COiT) = 




* 



CO© 



(A-11) 



* The bars indicate that the frequency coefficients are 
functions of ^ . 
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where 



P [(K) V y ) if ( r)*-//- y/ - 2o< (n*- /)*( *rr)2- y y V( o»- // '] ^ 



(A-12) 



- [(nVy) +®<y(n*-/f + '] 

w 

The frequency coefficients f^^ and may be Bluiplifled in Light of 

the thin ring assumption. For small values of eg *, these 
coefficients are approximated by 



CH 

yiHy 



(A-11) 



■y 



The initial values for the time functions must now be determined 
in order to solve for the constants end b^^. At time T * 0, the 
expresslnns for the displacements and velocities are simple Fourier 
series expansions in the slrtgle variable The Fourier coeffl* 

dents for these series may be evaluated to yield the foilowlitg 



Initial conditions for A^(T) and B (T) 

n It 

A (0) * 0 

ti 



n » 0, I, 2, 



1 • • 



*^?For the problem considered, o<^ 3.54 x 10 
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7T 






( 0 ) 




4 



^ Am / q\ - 2Vo(~0 

(Y)^-i)7T 



n 2jAj6j« • • 



= 0 



n = 3, 5, 7, 



B„(o)=0 




4 



( 0 ) = 0 



2, 3, 4, 



(A-14) 
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The time functions may now be determined from the above initial 



values. For the cases n = 0 and n 1, a simple substitution into 
the modified ring equations yields 

A, (T) = — ^ SlU COoT (A-15) 

7T Ode 



and 




The solution for the higher mode shapes, n > 1, requires further 
explanation. Recalling that the modal frequency COy^ was character- 
ized by four distinct roots, there are four separate solutions for 
and which satisfy the ring equations. The sums of these 
solutions are also solutions. Therefore, the general time functions 
for n > 1 are 



A -tfnOJoT _ - - 

Aw (T)= Q,r)i ® ^ ^ ■** 



I ^VtOJo/ 



(A-17) 



BniT)” b„, e' 



r i-^nCd^T r -^'^nO^cT 



■'h2 



r ^ I 






M 



Although eight constants appear in the above equations, only four 
are independent. As determined earlier in Equation (A-9) , the b 
constants may be replaced with expressions involving the "a constants. 
Applying Equations (A-13), these relationships are 
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n 


• Uv 


f) 









(A-18) 



— 


n 


• 


L - 7T 


n 


/ 


^K)3 




) 


Dr)4.“ ^r>4 







Substituting Equation (A-18) into Equation (A-17), the expressions 
for A (T) and B (T) contain four unknown constants to which four 
initial conditions may be imposed. 

To apply the initial conditions. Equations (A-14), the time 
functions must be expressed in terms of sines and cosines. Making 
the appropriate expansions for the exponentials and regrouping the 
a constants into new constants, a. Equations (A-17) become 



Ar> (T) = cos fr, Cd^T +• S//\/ 



>(A-19) 



+ cos l„0).r +nn^ZlN f, 

A solution may now be determined for the a constants. Applying the 
initial conditions, the cosine terms must vanish. Thus, 



- hfis = O 



(A-20) 



The coefficients for the sine terms are given by the following 
expressions 
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a 



r)4- 



2 Vo - ri^) 

2 V, 



2. y A, 6j •• 

(A-21) 



^|T)2 ^IT)4L - O 



rib- 



and hence 









-2V,(-i)”‘^*‘ n 
>rco,7T(io^-/)(f„2-5„')^„ 

irw.7T(ir)2-/)(f/-e^)^„ 

^K)2 ~ brt4 “ 



The time functions may now be written directly as 

A„(t)- 



n = 3 , 5 , 7 , . . 

(A-22) 



n = 2 , 4 , 6 ... 

(A-23) 



n = 3, 5, 7, 
(A-24) 











1 



2 vS2 






6,,(r)= 



2 V.±r^ 



)S0)^'n{ n^- /) (-^^ - 3^) 



SlN-PnCo^T +-§- S//V g«c0o7 

fyj 



even modes 2 

Aj^(r) = 6^fr)» <9 odd modes >3 



(A-25) 
(A-26) 

These expressions may be substituted into Equation (A-3) to obtain 
the middle surface displacements, W and V, as functions of time and 
ring position. As determined earlier in Chapter II, Equations (4) 
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and (5) , the stress-resultants M and N are given in terms of the 
displacements W and V by 






(A-27) 



Thus, the solutions for W, V, M, and N due to the radial cosine 
impulse become 



W ^ = SIN U.T - 



7T 






SIN ]Jl+^ Od^T Cos 4> 



Z(-0 



y>/2+i 












4- 



*^10 



SIN 



SIN -f^ CdoT 



COS n<p (A- 28 ) 



V 



CO^TT ^ 2 L 

V. 



4fJ7^ 



siN\j I ^ Y co^T Sin 4> 



oO 



+ 



Z(-I) 






n-z,4, 



-y) 



s/« 



siN-f„u)^r 



SIN n4> 



(A-29) 
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and 




A^tt g)o 

KVo 



oo 



= S/A/ co^r -t- 















S//\/ 4 



n 



-h 



Y)^ - /ftn 



S/A/ gn o^oT 



Cos n<^ 



(A-30) 



A4 = A/ 



AncOo 

bV. 



-S/A/ coj- V- 



7T 



S/W i// +k' WoT C<3S <#> 



oo 







2 (-iff*' 

"3«) 




S/N -fr)CO^T 



+ S//V c^n cOo T 



COS n<A - o< Mo 



(A-31) 
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APPENDIX B 



CALCULATED FREQUENCY COEFFICIENTS 
TABLE V 





FREQUENCY COEFFICIENTS 




n 


£ 

n 




2 


0.016 


2.236 


4 


0.088 


4.123 


6 


0.208 


6.082 


8 


0.376 


8.062 


10 


0.592 


10.050 


12 


0.857 


12.041 


14 


1.169 


14.036 


16 


1.530 


16.031 




TABLE VI 






FREQUENCY COEFFICIENTS 


( * 


n 


f 

n 


Sn 


2 


0.015 


1.732 


4 


0.085 


3.000 


6 


0.205 


4.359 


8 


0.373 


5.745 


10 


0.589 


7.141 


12 


0.854 


8.544 


14 


1. 167 


9.950 


16 


1.528 


11.358 
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TABLE VII 





FREQUENCY COEFFICIENTS 


( 




= 3) 


n 


f 

n 








2 


0.014 






1.528 


4 


0.083 






2.517 


6 


0.202 






3.606 


8 


0.370 






4.726 


10 


0.587 






5.859 


12 


0.851 






7.000 


14 


1.163 






8.145 


16 


1.524 






9.291 




TABLE VII 










FREQUENCY COEFFICIENTS 


( 




- 4) 


n 


f 

n 








2 


0.013 






1.414 


4 


0.081 






2.236 


6 


0.199 






3.162 


8 


0.367 






4. 123 


10 


0.583 






5.099 


12 


0.848 






6.082 


14 


1.161 






7.071 


16 


1.521 






8.062 
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